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SYMSIS 

HYBRID SIMUMTIOH OF THE SUPERHEATERS OP 
THE 110 W BOILER UNIT OF PANKI 
THERMAL POWER STATIOH 

A Thesis Submiuted 

111 Partial Fulfiliaent of the Requirements 
for the Degree of 
MASTER OF TEOHHOLOGY 

By 

K. Balagopal 

The non-linear mathematical model of the superheaters 
of 110 MW Panki Thermal Power Station is simulated on hybrid 
computer system. The responses are compared with digital 
simulation responses. It is seen that the maximum error 
between the results is of the order of 1 percent, A linearized 
form for the non-linear model has also been obtained and simu- 
lated, The results compare well with the non-linear model. 

The advantage of using hybrid simulation over digital compu- 
tation can be had from the fact that the computation time 
required in the hybrid mode for the non-linear model was only 
5 minutes (at 1/6 slowed down scale) whereas a 4th order 
Runge-Eutta digital simulation on the sane computer (IBM 1800) 
in digital mode took two hours. 
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IITRODUCTIOH 

1.1 HYBRID COMPUTERS - AH IHTRODUCTIOH 

In the pursuit of accuracy, speed and economy in the 
simulation of mathematical models of physical systems, ana- 
lysts are continuously searching for more and more sophisti- 
cated mathematical techniques and computer systems. When 
the complexity of the mathematical model increases, their 
simulation on pure analog or pure digital computers are 
frequently cumbersome. The main reason for this is that 
some of the operations in mathematical model are inherently 
best handled either on the digital or on the aiaalog computer. 

The chief distinction between analog and digital compu- 
ters lies in the manner in which the dependent variables are 
handled within the computer. In analog machines, dependent 
variables (though not necessarily the independent variables) 
appear everywhere in continuous form and may be recorded with 
as many significant figures as the quality of the circuit 
component permits. In digital computers, all variables 
appear in discrete form, and the accuracy of the data which 
are manipulated and recorded as the output of the digital 
computation depends on the number of digits carried through- 
out the solution and is directly related' to the capacity 
or size of the memory registers. 
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Some advantages of analog computers are, 

(i) the dependent variables are in continuous form^ 

(ii) the operation is parallel, with all computational, 
elements operating simultaneously and hence high-speed 

or real-time operation, (speed depends upon the bandwidth 
characteristics of the elements, not by the complexity 
of the problem)^ 

Some of the draw-backs of analog computers are, 

(l) lesser accuracy, accuracy depends upon the quality of the 
computer components and rarely better than 0.01 percent, 
(il) very limited ability to make logical decisions, store 
nxamerical data and provide extended time delays. 

The points in favour of digital computers are , 

(i) facility for memorizing numerical and non-numerical 
data indefinitely , 

(ii) facility to perform logical operations and decisions, 

(iii) accuracy does not depend much on the quality of the 

system components, but primarily depends on the number 
of bits contained in memory registers and by the spe- 
cific numerical techniques selected for a particular 
problem. 

The main drawback is that the computer time is relati- 
vely long and depends on the complexity of the problem 
(serial operation) involvirg the time sharing of all operations 
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and memory units, only one or a limited number of operations 
being carried out at one time. 

The Hybrid computer technique is an effort to combine 
in one computer system the advantages of both the analog 
and digital computers. The chief motivations behind hybridi- 
zation are to combine the speed of the analog computers 
with the accuracy of the digital computers, to increase fle- 
xibility of analog simulation by using digital memory and 
logic capabilities, to increase digital computation speed 
by utilizing parallel analog subroutines and to handle data 
in systems which are partially continuous and partially dis- 
crete. 

In a system of differential equations, if the frequency 
of solution of one variable is considerably different from 
the frequency of solution of the other variables, then it is 
very difficult to solve the equations efficiently with any 
numerical method. The same step size will not be suitable for 
all the variables. Pure analog simulation may also be 
difficult depending upon the complexity of terms. In such 
cases hybrid method can be effectively employed. 

Hybrid computers are very effective for multivariable 
function generation. Hardware of costly non-linear elements 
are required to generate functions of two or more variables 
in an analog computer. Digital computers cannot be employed 
for generating continuous functions. 
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Hybrid computers are more efficient for solving partial 
differential equations compared to digital or analog computers. 

Some of the major fields in which hybrid computation is 
being used are simulation of aircraft vehicles, missiles 
and space vehicles [9J. Hybrid simulation is used for the 
design and safety study of Huclear power plants, 

A hybrid computer system in its simplest form consists 
of an analog computer, a digital computer and an interface 
unit between the two. The interface should consist of a 
multiplexer, analog-to~digital converter (ADO), digital-to- 
analog converter (DAG) and a demultiplexer for communication 
between the two computers, a timing and control unit which 
synchronizes and controls the operation of the hybrid loop 
(the timing and control unit is mostly part of the digital 
computer). The block diagram of a hybrid computer system is 
shown in Fig, 1,1 
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Fig, 1.1 Block diagram of a Hybrid Computer 
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1.2 HYBRID SYSTEM AT I.I.T. K. 

Bor the hybrid system, the digital computer available 
was IBM 1800 (Data acquisition and control system) and the 
analog computer x^as AC 20 of ECU make. A process computer 
should have the analog input - output terminals, process 
interrupt facilities, and interval timers IBM 1800 has these 
facilities and can be used as the digital part for a hybrid 
computer system. 

The specifications of IBM 1800 and AC 20 are given in 
Appendix A. 

Choudhary [3] and Sastry [14] of Electrical Engg. 
Department and Eaghavan [10] of Computer Science programme 
have previously worked on the present facility. 

1.3 OBJECTIVES 

The main objective of the present work is to -employ hybrid 
computer system for the simulation of thermal power plant, 
specifically the 110 Wil unit of Panki Thermal 
Power Plant, ^ince the hybrid system at I.I.T.K. is very 
small, large systems cannot be simulated on this. So, to 
start with, a small sub-system - the syperheaters of Panki 
Thermal Power Plant was selected for the present work since 
the number of dynamic equations are well within the limits 
of the present facility. A mathematical model of a Boiler 
unit is given in the IBM Report [5]. Sreenivasan [16] 
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ier eloped a mathematical model for the Boiler unit of Panki 
Thermal Power Plant similar to this model. The mathematical 
model of the superheaters was selected for the present work. 
The model of the superheaters consists of seven differential 
equations only which are nonlinear and coupled in nature. 

In the present work, the nonlinear model has been lineari- 
sed around an operating point considering small perturbations. 
The hybrid simulation results for both linear and nonlinear 
models have also been compared. 

One another objective is to find out how much accuracy 
the present system will give in simulating real systems. It 
is thus necessary to compare the hybrid simulation results 
with the actual resxxlts, here the pure digital simulation 
results using fourth order ibinge-Kutta method. 

De Boer and Dehes [4J simulated the steam generating 
system of a IMPBR with control loops in a hybrid computer 
system. They used applied Dynamics - 4 /IBM 1800 hybrid com- 
puter facility at the Delft University of Technology for their 
work. The analog computer was used for the continuous inte- 
gration and for the calculation of some of the parameters. 
Digital computer was used to calculate the coefficients used 
in the analog computer, to play back data from the previous 
time step to the analog computer, to calculate some of the 
parameters and to store data. 
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Paul landaner [9] reported the simulation of the 
Pickering Nuclear Gtenerating Station control to supplement 
the design work of individual loops and to review the per- 
formance of the plant control system before the plant went 
upto power, 

1.4 fODAI'S HYBRID OOfrlPUTER [13] 

On sharp contrast to the relatively crude combinations 

of largely manually set analog computers alongwith second 

generation digital computers, today's hybrid computers are 

Tn uttiprocA^^ofS 

fully automatic, preprogrammed with problem 

Set up accomplished in seconds and new parameters inserted 
and results obtained in micro or muliiseconds. Thus systems 
are designed around advanced, hybridized analog and digital 
circuitry, together with software that can automatically 
provide the optimal combination of analog and digital solu- 
tions, operator-interactive high-level application of pro- 
gramming languages, and disk oriented foreground/background 
real-time operating systems. The average hybrid computer 
system has a speed equivalent to more than 100 million 
digital operations per second, giving a speed advantage of 
20 to 150 times that for digital-only solutions for dynamic 


models. 
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aEiPISR 2 
SUPERHEATER MODEL 


2.1 INTRODUCTION 

For tile reasons mentioned in Chapter 1, the superheater 
part of the Boiler unit of the Panki Thermal Power Plant was 
taken for the present work. The boiler is a single drum, 
radiant, water tube, dry bottom circulation, coal fired 
boiler with superheaters and reheaters. Some important 
specifications of the boiler unit are given below [10] . 

Bfenufacturers ; Bharat Heavy Electricals Ltd., India 
Full Load Specifications : 

Generating Capacity 110 MW 
Steam flow rating 375 Tons/hr 

live Steam temperature 540° 0 
(At final superheater outlet) 

Pressure 139 Eg/cn^ 

(At fina.'. superheater outlet) 

Drum Pressure 150 Kg/cm^ 

Fuel Pulverised coal 

Calorific value of the fuel 32 percent by weight 

The general layout and description of the plant is 
given in reference [10]. 



9 


Steam is collected through, the 5 epa- 2 :*ators from the top 
position of the drum (Fig. 2.1). This steam at saturated 
conditions is led to ceiling superheater through headers. 

The steamf coming out of the ceiling superheater is passed 
on in turn to convection superheater., platen superheater, 
and final superheater respectively. This superheated steam 
is supplied to steam tTU'hines through throttle valves. The 
temperature of the superheated steam is controlled hy attem- 
peration at two places, i.e. between convection and platen 
superheaters and final superheater. The physical location 
of the various superheaters in relation to the furnace is 
as follows. Ceiling superheater comprises of the exposed 
surface of the furnace, platen superheater is situated at 
the top of the combustion chamber, and the final super- 
heater at the exit of combustion chamber. Ceiling, platen 
and final superheaters receive heat mostly by radiation 
whereas convection superheater sit-uated in the convection 
zone receives heat mostly by convection. 

The superheaters are represented in three sections. 

The ceiling and convection superheaters are clubbed together 
and is termed as the primary superheater. The second 
section is assumed to be a unit which effectively simulates 
the flow dynamics of the superheater lines. This is imagined 
as a connecting link between primary and secondary super- 
heaters. The secondary superheater consists of platen and 



Primary super Connecting Secondary super 

heater link heater 



Fin 2 I Model reDresentation 
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final superheaters. It receives steam flow from the pri- 
mary superheater through the linking line and discharges 
steam to the turbine through the throttle valves. The 
temperature of the superheated steam is controlled by 
attemperation carried out in iDrimar^y and secondary 
superheaters. 

2.2 PRIMRY SUPERHEATER 

The mass and energy balance equations of primary 
superheater are given by [Section A to B of Pig. 2,1 j. 


(m^i - W; + ) 


DSl- 


^^Sl ^Sl^ 


dhgi dqgj_ JPg^ d ^32. 1 

-dF" = ^ XrFTITjJ 

si 


( 2 . 1 ) 

( 2 . 2 ) 


where 

m^^ = evaporation rate or steam flow rate leaving the 

drum in kg/sec, 

¥' = mass flow rate of steam leaving the primary super 

0 

heater in kg/sec, 

Psi = density of steam in the primary superheater in 

kg/np , 

3 

V = volume of the primary superheater in m-'^, 

B JL 

¥-noT = mass flow rate of o, t temp era t ion spray added in 
Dbl 

kg/sec, 

J = conversion factor (work to heat energy) 
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= net heat content in the prinary superheater per 
kg of steam in kcal/kg and is given by 


_ ,-r a. V, 

dt 1- si ®ri wos 


iflf^h .j_ + ¥pg_^ hjjg] 


P sl^sl 


(2.5) 


where 



h 


wos 


h. 


DS 


h 


si 


si 


= external heat added to the steam (primary superheater) 
in kcal/sec. , 

= enthalpy of saturated steam entering primary superheat ea 
at drum pressure in kcal/kg, 

= enthalpy of the attemperations spray in kcal/kg, 

= enthalpy of steam in primary superheater in kcal/kg, 

= pressure in the primary superheater in kg/m® and is 
given by. 


^S1 = ^S1 “ ^2^ ^ SI 

where and K 2 Oallendar's proportional constants, 

2.3 SECONDARY SUPERHEATER 

The mass and energy balance equations of secondary 
superheater are given by [Section C to D of Eig, 2,1J, 



W"t + W- 


DS2^/’^S2 


d 

dt 



(2.5) 


*^^82 _ .^’^ ^32 
dt ~ L dt 


JPS2 

PS2 


^ ^32-, 1 

“dt -i r K^j) 


( 2 . 6 ) 
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where 


¥. 


DS2 


g 2 = density of steam in secondary superheater in 

mass flow rate of attemperation spray in kg/sec, 

'Z 

"^32 ~ volume of the secondary superheater in m"^, 

= mass flow rate of steam going out of superheater 
through a throttle valve in kg/sec, and is 
calculated using the relation, 

(2.7) 


- GA (Pg2 Psp) 


where 


0 = throttle valve coefficient, 

2 

A = area of the throttle valve in m , 

2 

Ps 2 = pressure in the secondary superheater in kg/m 
and is calculated using the relation. 


^S2 = (% - 22)Pb2 C2-8) 

where 

hg 2 = enthalpy of superheated steam in the secondary 
superheater in kcal/kg, 

q„o = iiet heat content in the secondary superheater in 
kcal/kg and is given by. 


dq 


s2 


dt 


= [Q, 


's2 




%S2 *^"^3 




where 



14 


Qg 2 “ ©;^ernal heat added to the steam in kcal/sec (secondary 
superheater) . 

2.4 CONNEailNG LISK 


This section (section B to 0 of Fig. 2.1) is assijined to 
be a unit effectively simulating the flow dynamics in the 
superheater lines. It can be thought of as a line connecting 
the primary and secondary superheaters which has no net heat 
loss or gain. The pressure drop across this \mit determines 
the flow between the primary and secondary superheaters and 
this flow is modified by first order dynamics to incorporate 
an approximate representation of flow dynamics and is given 


d¥' 

jf- = S (Ssi - Psa) - ''i 


( 2 . 10 ) 


where 


T 


GL 



= time delay constant for the connecting link in 
seconds, 

= a constant determined at the time of initializatioK 


The time constant T^^j^ introduced will* account for the 
dynamic lag in flow that would be experienced with changes 
in pressure and the transport delay that is associated with 
the superheater line lengths. 


2.5 SUPEEHmTBR TUBBS 

The energy balance equations for the primary and 
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secondary superheater tubes are g-iven by 


(QgMl - 


“ dt^^SICL^ 

(2.11) 


«S2)/°S2«2 

11 

(2.12) 


where 


Cgi, 0^2 = specific hO-.t of primary and secondary super- 
heater metal respectively, (kcal/kg'^C) 

Ml, M2 = lotal weight of primarysupsuperh eater metal, 

secondary superheater metal respectively, (kgs) 
'^SM 1 ^%M 2 ~ Heat transfer rate to primary and secondary 

superheater metal from gases respectively, 
(kcal/sec) . 


Ihe rates of heat transfer from the tubes to the steam are 
calculated using empirical relations. Since the flow is 
turbulant, the relations used are 

“^sl %s ^^Sm " ^sl^ (2.13) 

^^82 ^^SM2 “ ^s2^ (2.14) 

where 

Gpg,Gj^g = Coefficient of heat transfer from superheater 
metal to steam in primary and secondary 
superheater respectively. 
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“ lurbulant heat transfer index, 

'^SM1*^SM2 ~ IsEperature of the prinary and secondary super- 
heater metal respectively, (°C), 

^sl»^s2 “ Temperature of s’- earn in -hie ]prinary and secondary 

superheater respectively, (°C) and are calculated 
using the relations , 

^sl ~ %l^^sl “ ^sTl^ ^Rl^ Psl ~ P sTl^ ^sTl (2.15) 

^s2 ~ %2^^s2 “ ^sT2^ ^R2^ Ps2 “P sT2^ ^sT2 (2.16) 

Tgji, Tg ^2 = Steady state operating temperature of steam in 

the primary and secondary superheater respectively, 

^sTl’ ^sTl' ^sT2 P sT2 ^©P^^esents the values h^^, hg 2 

and Pg 2 respectively at steady state operating temperatures. 
Now, the differential equations can 'be rearranged with seven 
isRdependent variables as follows ; 

It ( I’si) = ('“ri - «dsi)A 31 (2-in 


dt 


= t[ %1 * “ri \oB - '^o ^sl ^Dsl 

+ I 1 

Pg2 >1* ■' (1 - Kp) 

If ( fsa) = - “l + %S2 ]i4 



Psl "^Sl 


(2.18) 


(2.19) 


dt^^s2^ ~ t 


Qs2 


+ 

4- 




jp 


s2 


Pg2 

®2 


^Ds2 ^ds ~ ^1 ^s2]pg2^g2 

It ^‘^8 2^3 (i" r 


(2.20) 
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K. 


fW’ ) = —£• 


OL 


d?' ^^sMl “ *^31^ Cc ]_Mi 


d _ 
dt ''■^sM2 


~ (^PlM9 ~ ^.ejo) 


1 

'SM2 “ '^8 2^ *^3 2^2 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


2.6 LINME MODE! 

2,6.1 Linearization [15] 

In many investigations of nonlinear systems, a solu- 
tion for a given set of initial conditions and excitation 
has been determined (either analytically or with the aid of 
a digital computer) and it is required to determine the 
behaviour of the system in the neighbourhood of this solu- 
tion. Eor small perturbations from such a solution, the 
a^nalysis reduces to that of a linear system. 

Consider a nonlinear system with initial state 
X(t^) and input U(t) in the form, 

X = f(X,L,t) (2.24) 

The Solution of this equation is X(t). 

Let the initial condition is changed to 2(tQ) + x(t^) 
and the input to u(t) + U^(t), where x and are both 
very small. 
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Since the resulting state will satisfy the differential 
eq.ua t ions y we can write 

i + z = f(X + z, U + U^, t) (2.25) 

Bzpanding the right hand side of (2.25) in Taylor's series, 

• • -P Ci -P 

Z + z = f(Z,U,t) + z + ^ + higher - order terins(2, 26) 

(i.e.) z = 1^ z + 1^ + higher - order teras (2.27) 

If the higher-order terms are small, they may be 
neglected in obtaining an approzimate solution to the 
problem and the resulting perturbation equation is seen 
to be linear, 

2,6,2 linear Mod el 

Using the above method, all the non-linear equations 
are linearized in the present work. The state and input 
vectors are as follows ; 

~ t ^sl ^sl ^ s2 ^s2 ^o ^sMl ^sM2^ (2.28) 


= CA ¥j,g2] 


(2.29) 


As an ezample, the linearization of eqn, (2.19) is given 
below. 


d 

dt 


(W«) = 


K 

JSL. 

^01 


(P 


Si 


I’s2) 


W 


0 



p 



19 


where 


^sl “ " ^2^ ^sl 

= (^1 1^S2 - ^2) *’ "2 


After substitution eqn, (2.19) can be written as 


VA .£x.q xa. 

ar = [»1 ^2 - ^z') ^4 - - 

linearizing, 

dXc - K 

aT = I?. ^2 ^10 + (^1 ^20 - ^1 - 


T 


GI 


•01 


[% ^4 ^50 ■" ^40 - ^ 2 ) ^3^] “ T 


fi 


01 


where X^q> ^20 and ZjQ represent the steady-state 
of Xj^, X2 and X-^ respectively. 

This can be written in the form, 

= A(5,l) X 2 _ + A(5,2) ^2 + A(5,3) z^ + A(5,4) z^ 

A(5,5) Z 5 + A(5,6)x 5 + A(5,7) z^ + B(5,l) 

+ B(5,2) U3_(2) -h B(5,3) %(3) 

where the coefficients are given by 

A(5#l) = (% X 20 “ ^2^ 

K 

A(5,2)=-ffi^ (E, S,o) 


01 


( 2 . 30 ) 


(2.51) 

values 


( 1 ) 

(2.32) 
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A(5,3) = - ^ (El - Kg) 

Oil 

E 

A(5,4) = - ^ (% X^q) 

GL 

A(5,5) - -1/^cL 


A(5,6) = 0 
A(5,7) = 0 

B(5,l) = B(5,2) = B(5,5) = 0 

Similarly all the equations are linearized and the model 
obtained is in the form, 


X = [A] X + [B] TJ 

The operating point taken for the linearization are given 


below. 


i = tA WbSI ''dS2] = 5.472] 


- 


11A327 

hq 

^1 


663.786 

^S2 


40.1474 

hg2 


819.112 

¥' 

— 

107.354 

0 



^SMl 


380.099 



568.09 


The matrix A and matrix B obtained after 


linearization are 


shown below t 
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CHAPTER 3 

SOFTWARE AED SET UP FOR HYBRID OO^IPUTATION 

3.1 IETaODUCTION 

Since an overall idea about tlie IBM 1800 TSX system is 
necessary for the understanding of the softt'/are used, it will 
be discussed briefly first. Although Ragimvan [10] in his 
work had developed the software required for the hybrid 
system, a new set of subroutines are developed in this case. 

This was necessary because one of the main tasks was to 
reduce the digital computer time as far as possible. A lot 
of variables and subroutines used in Raghavan' s software 
package were unnecessary for the present problem, 

3.2 IBM 1800 TSX SYSTEM 

Time inhering Executive (TSX) and Multiprogramming 
Executive (MPX) are supervisory programs supplied by IBM to 
operate IBM 1800 computer. MPX is an improved and sophisti- : 
cated version of TSX. Since to use MPX efficiently the compute] 
should have at least 24R memory, the IBM 1800 at I.I.T.K. 
uses TSX. TSX is 4 FORTRAN oriented disk - resident opera- 
ting system which permits the user to make optimum use of 
an IBM 1800 Data Acquisition and Control System (DAGS) for 
its primary purpose, the control of processes and similar 
environments, as well as providing the user with an effective 
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off-line monitor system for data processing and scientific 
computation. 

A 'process program' is one that continuously monitors 
a control process, A process program uses process subroutines 
through which the process input - output devices are operated 
like digital and analog input - outpu.t devices, interrupt 
device, clocks etc. For the hybrid computation, the process 
capabilities of IBM 1800 are used to give analog signals, to 
read analog inputs and to control the modes of the analog 
computer. Process programs are classified into three t 

i) Interrupt Program 

This is a process program which can be executed by 
creating an interrupt. An interrupt can be created by one of 
the following methods a) pressing the 'EBD REQ' key of the 
l/O typewriters, b) closing of a process interrupt contact 
point (these points are displayed on a panel on the back- 
side of the IBM 1800). 

ii) Mainline (Queue) Program 

This is a process program which is executed by queuing 
it. A mainline program can be placed in the queue by any 
process program. 

iii) Combination Program 

This is a process program which can be executed either 
by queuing it or by creating an interrupt. 
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3.3 MODE COITROI OF AMIOG COMPUTER 

In a hybrid set up the three main inodes of the analog 
computer, RESET, COMPUTE and HOID are to be controlled by 
the digital computer. Thre subrouiiines written serve this 
purpose. This is actually attained by using the digital 
output ~ ECO (Electronic Contact Operate) points. The relay 
positions of integrators for mode control [1] and the circuit 
designed by Choudhary [3] to attain it using the IBM 1800 
are given in Appendix B, 

3.4 ANALOG READ - WRITE 

Two subroutines (other than the system subroutines) 
are necessary for analog read and analog write because of 
the disparities in voltages between analog and digital com- 
puters. The analog computer voltage varies from -10 V to 
+ 10 V. The analog output voltage of the digital computer 
ranges from 0 to +10 V for the outpuc points 00 and 01. 
Remaining output points give output in the range 0 to 
+ 5 Volfs. The analog input to the digital computer can 
vary from -5 V to +5V (This is when solid state multiplexer 
is used. For relay multiplexer it is 0 to + 5 volts,- only 
one sign is applicable at a time). Level shifting circuits 
are necessary to compensate these disparities. The 
circuits used are shown in Fig. 3,1. 
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3.5 HIBHID ROUTIBBS 

Program RAGVN 

This is an interrupt coreload program. When the 
interrupt button on the interface unit is depressed, this 
routine will be transf c:'“red to cor' and executed.. This 
program queues BAIGO. . 

Program BALGO 

This puts the analog computer in reset node and 
initialises all output voltages to zero. This queues the 
hybrid program 'HYBR' . 

Program RAGS 

This is a mainline coreload program. This is a 
dummy routine and is necessary to delete HYBR. (HYBR is 
called in BAIGO and it can not be deleted unless its call 
in BAIGO is substituted by another mainline coreload pro- 
gram name. RAGS serves this purpose). 

Subroutine RESET 

This puts the analog computer in reset mode. 

Subroutine COMPT 

This puts the analog computer in compute mode. A timer 
is called in this routine immediately after the analog com- 
puter is put in compute mode. The calling statement of this 
routine is CALL COMPT (TYM). ' TY14' specifies the compute 
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time of the hybrid system in seconds. The routine ‘TBCS’ 
will be called after every 25 seconds and a variable is 
equated to 1, the value of which was initially zero, when 
the compute time elapses. So, to check whether the compute 
is over or not, the value of this variable is checked in the 
hybrid program frequently. If this value is 1, the analog 
computer will be put in reset mode, otherwise, the computa- 
tion will be continued. 

Subroutine EDA I 

This reads sequentially 7 analog input points to the 
digital computer and converts them to the corresponding 
analog voltages starting from the address point K and puts 
the values in the array X(l) to X(7). The calling statement 
is CALL EDAI (X,E). 

Subroutine ¥RAO 

This first makes the voltage 'VAL'tFB' to be written in 
the digital computer output range and then writes the value 
at the particular address point ' JPT’ . The calling statement 
is CALI ¥RAO (VALUE, JPT). 

Subroutine TIME 

This routine is called every 25 seconds. ¥hen the 
compute time is over, the variable 'BALU* in this routine 
is equated to 1. 
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Subroutine TIICBF 

Tills is to fool the corapiler of IBM 1800 because 
it does not permit the rame of a subroutine within itself. 
Routine -TIME' is to be cal_ed ewgry 25 seconds and this 
is not possible within 'Til®' itself. So Til®! is called in 
TIME and this calls TIME every 25 seconds. 

Subroutine SST 

This subroutine is called in 'HYBR' freq.uently to check 
whether the compute time is over or not. 

5.6 SEQUEEOE OE EXECUTION 

When the interrupt button on the interface unit is 
depressed, EAGfVN will be executed. This queues BALGO. \Ihen 
a program is queued the programs existing before the begi- 
nning of the execution of the queued program are overwritten 
in memory and this program is executed whereas a called 
program always forms part of a coreload. BAKiO queues the 
hybrid program 'HYBR* . The flow chart of the hybrid program 
is given in Pig. 3.2. 

3.7 PROBLEM PATCH-UP OH HYBRID SYSTEM 

Out of the seven differential equations, the first 
one (eqn, 2.17) is solved o©mpletely on the analog computer. 
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The operational amplifiers and potentiometers of AC 20 are 
used for the problem patch~up. The potentiometers and 
operational amplifiers of TR 20 are used for the level shifting 
circuits. In general, the a:ialog part is u^ed for the conti- 
nuous integration and for the level shifting circuits. The 
digital port controls the modes of the computer, reads the 
analog voltages, calculates the algebraic equations and gives 
back the signals to the analog computer and stores the inter- 
mediate values. The complete set up is shoun in Rig. 3.3. 

3.8 SGALIhG THE EQUATIONS 

3.8.1 Time Scaling 

Time scaling is used in analog computers to slow down 
or Speed up the solution. The idea is that the maximum fre- 
quency of oscillation of the solution of a differential 
equation should be within the bandwidth of the analog com- 
puter components. 

In hybrid computers time scaling is used to slow down 
the analog computer so as to have the required inmaber of 
samplings in one cycle of the highest fresiuency variable. 

One of the most serious errors in hybrid system is the 
time-delay error. This is the effect of digital execution 
time in a hybrid computer. The analog to digital conversion 
of the voltages, its processing and converting back to analog 




.3.3 Problem pa 
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voltages will take a finite aaomat of time. During this 
time the analog computer would have continued to compute. 

So, if the amount of computation on the digital computer is 
increased in order to raise accuracy, execution time is 
prolonged and so, increased computing error may resiilt. 
According to Tr.t;it% |[l8]j ^ good rule of thumb, a sampling 

rate of about 50 to 100 per cycle of the highest frequency 
variable will keep the time-delay error within 0.1 percent. 

Prom the linear model of the present system, it oin be seen 
that the frequency of the highest frequency variable is , .. 
3,65 rad/sec, (i.e,) 0.581 cycles/sec. In the real time 
simulation the data were transferred 120 times in 50 seconds, 
that is 2,4 times in a second. This shows that in real time 
simulation, the time-delay error will be considerable for 
this system. The analog computer was slowed down using 
= xJS for the nonlinear model and t = t/ 4.5 for the 
linear model [11] » 

3,7,2 Amplitude Scaling 

The purposes of amplitude scaling are, i) the amplifiers 
should not get overloaded, ii) the full dynamic range, namely 
+ V (the range over which the amplifiers are designed to 

work — reference voltage) of all amplifiers should he utilized* 

In the present work, to utilize the full range of all the 
operational amplifiers ^ the origin of tno state variables 
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was shifted to a new value, lliis is done by using 

Z = X - G (3.1) 

where Z represents the new state after shifting the original 
state X ■:o the origin G. All the equations were reformulated 
using this scheme. The values of the origin G for various 
percentage openings of the throttle valve area are given in 
Table 3.1. 

After shifting the origin the equations were amplitude 
scaled using the normalised variable technique [11]. 

A listing of the subroutines and programs used for the 
present hybrid simulation is given in Appendix 0. 



Table 


35 


0 





4> 

-H > 

03 



0 

0 

0 O 


n 

m H 

•H 

U 

d 4^ H 

0 

0 -P 

H 

ft 

M O 

c5 ^ 


o U 

0 tJ 

O 

^ xi 

JLi O 

H 

•H -P 

g3 si 


I 

I 

I 


I 


^ CM CO 

} i i 


CM LTN 
4- I I 


i 


I 


® 

j> 

H 

d 

-P -H f> 

r-j 

0 

[ o 

1 

1 

1 

£> 


® ® © 

a S 

1 

CM 

CD 

o ca H 

o 


{A 

00 

d 

Jit 

1 


C-- 

© © +=> 

d 


• 

* 

ft f4 O 

g5 2 

1 

H 

!> 

O fs 

0 d 


00 

CO 

O © 
iHf -P 

f-i 'H 

0$ H 

i 


CO 


t>- 

CM 

rH 

EA 

cr> 

cn 

O 

C3^ 

H 

H 

!A 

o 

O 

« 

• 

• 

♦ 

• 

CM 

tr- 


CM 

tA 


CM 

O 

CO 

CD 


CO 

H 

lA 

CO 


H 


0 


1 



l> 





H 

ft 

1 



^ c5 

0 




-P -H > 

ft 

r 



a 

1 O 




0 0 0 

^ a 

1 



O 0 H 

o 


CM 

cn 

«H Cd “P 

PJ u 

i 

A 

00 

0 0 -P 

c5 



{>- 

ft o 


T 

* 

ii 

o u 

0 d 


crs 

LA 

o 

ft 

1 

CO 

in 

CM -H 4^ 

03 ft 



CO 


ft" 

CM 

ft 

A 

(T\ 

<y> 

o 

cr\ 

ft 

ft 

EA 

o 

o 

« 

f 

• 

• 

» 

CO 

tA 

m 

CO 

00 

bA 

O 

ft 

A- 

in 


00 

ft 

EA 

in 


i 

} 0 

d d 
-P -H i> 

1 M 

10 0 0 

ft 

0 

ft 

I o 

5^ a 

I 

1 

r 

i 

t 

CM 

CO 

ft- 

i>- 

CM 

ft* 

cn 

o^ 


i o 0 ft 

o 


tA 

CO 

ft- 

ft 

^A 

o 

CA 

^ 03 -P 

^ u 

j 

ft- 

o- 

ft 

ft 

lA 

o 

O 

i 0 0 -p 

d 


• 

9 

« 

« 

• 

« 

• 

ft p 

d 0 

i 

EA 

<n 

CO 

ft 

O 

03 

A 

o u 

0 a 



in 

A 

ft 

ft 

!>* 

CO 

o f::? d! 

ftft-P 

U ft 
g 3 ft 

i 

1 

t 


CO 


CO 

ft 

A 

in 


H H CM H CM 

m m 0 CM -o S S 
CL ;ci ^ 0 :b= CO cq 

^ ^ Si 


1 


i 


I 

I 

1 



36 


GHAPTBR 4 

RESULTS AND DISCUSSION 
4.1 RESULTS OF HYBRID AND DIGITAL SIICJL^LTION 

Hybrid simuJ-ation of a themal power plant sub- 
system lias been done in this work. The sub-system consists 
of the primary and secondary suiieifieat ers and their linking 
component. The highly nonlinear mathematical model has been 
linearized around a normal operating point of the plant. 
Simulation is done with both the non-linear and linear model. 

The digital simulation results were obtained using the 
4th order Runge-Eutta method. In hybrid simulation, the 
analog computer was slowed down using t = t/6 for the non- 
linear model. So, for 50 seconds of plant operation, the 
time taken in hybrid simulation is 500 seconds or 5 minutes. 
It is to be noted that the digital simulation using Runge- 
Zutta method took two hours on IBM 1800 to get 50 seconds 
of plant response (It took 'nly 6 minutes on IBM 7044). 

Using the time scale t = t/S data were transferred 
722 times between analog and digital computer within 50 
seconds of real time (the plant operation time), that is, a 
sampling rate of 14.4 times/sec. or 24.6 times/cycle of the 
mazimum frequency variable. This sampling rate is suffici- 
ently high and the time delay errors can be ignored. 



37 


In tlie linear models since tlie complexity of terms in 
the digital computer were less, a higher sampling rate X'jas 
used. We obtained a sampling rate of 14.16 tirnes/sec. with 
t = t/ 4.3. That is the hybrid simulation of linear model 
took only 225 seconds for the 50 seconds plant response. 

This clearly shows that if we can reduce the calculations 
done in digital computer, the solution tine can also be 
reduced considerably. 

The responses obtained for the various perturbations 
were satisfactory. A step increase in throttle valve area 
(Pigs. 4.1, 4.2 and 4.4) means an increase in steam flow- 
rate, w^. For the constant attemperation sprays, the 
step increase in throttle valve area reduces the densities 

and Pg2» Similarly since the heat input is constant, 
the enthalpies hg^ and hg2 and temperatures Tg^^^^ and Tg5,^j2 
decrease. All these stabilise ai'ound a new operating point 
after some time. 

Similarly a reduction in throttle valve area, with the 

other inputs as constants reduces the steam flow w' and 

o 

increases all the parameters, Pg^^ ^Sl’ ^ S2* ^S2* ^SMl 

^SM2 ^*3). 

The maximum error between the results obtained from 
hybrid simulation and that from the digital simulation is 
seen to be only of the order of 1 percent (Figs. 4.1 to 4.4). 
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4.2 LIMITATIONS OF THE PfiBSEIvT HYBRID FACILITY AlCD 

SUGGESTIONS FOR FURTHER WORK 

In the present hybrid facility only a uaxinun of ten 
integrators are available fo:-. use. large sjsten nodels 
having more than ten differential eG'Jations co,nnot be simulated 
in this Set-up. If the analog computer capacity can be in- 
creased, larger systems can be simulated. Presently only 
one AC 20 and three TR 20 are available and the AO 20 is 
coupled to IBM 1800. If one TR 20 is coupled to the IBM 1800 
the other two TR 20 can be slaved along with this so that a 
maximum of 24 integrators will be available for simulation. 

One other limitation is that the operational amplifiers 
of the analog computer are subject to drifts if operated for 
a longer time. In simulating larger systems, the digital 
computation time is increased but the sane sampling rate is 
to be maintained to avoid time-delay erroi-s. In order to 
match the increased digital execution time, the analog com- 
puter has to be slowed down by time-scaling- Thus the analog 
computer will have to be in compute mode for a longer time 
interval where by the results nay be affected due to the 
drifts as mentioned earlier. 

The IBM 1800 at I.I.T.K. can give only 6 analog outputs. 
For simulating larger systems we need more than six analog 
outputs . 
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If we can use separate digital-to-analog converter in 
conijmction with. IBM 1800, the digital output terminals can 
be used to give analog signals thereby enhancing the number 
of analog outputs. 

During the simulation of larger systems, digital or 
analog compensating techniques [17] can be used to reduce the 
time-delay errors. Uo such compensating method was used in 
the present study because a suitable choice of time - scaling 
was possible and which gave satisfactory results. 
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APPEITDIZ 


A 


A-1 IBM 1800 Specif i cat ions [6] 

1. Word length : 18 bits (16 data bits, 1 bit for parity 

check, 1 bit for storage protection) 

2. Main Memory : 16 2 words (expandable to 64K words) 

3. Memory Oycle Time : 4 psec. 

4. Instruction execution times : 9 psec. for addition, 30 psec, 
for multiplication, 85 psec. for division etc. for integer 
arithmetic. 

5. Analog Input 

i) Solid State Multiplexer - Switching rate is 10,000 
points per sec. 64 analog inputs can be connected. 
Permissible voltage range is from -5^ to +5^. 

ii) Relay Multiplexer - Switching rate is 100 points/sec, 

15 analog inputs can be connected. Permissible 
voltage range is from— .5 to +5 V. 

Digital input 

12 sets (16 bits each) of digital input terminals. 

Out of these 12 sets, two sets are of contact sense type 
and the rest are of voltage sense type. 

7. Analog output 

Six analog output terminal points (addresses from 00 
to 05) 00 and 01 give voltages from 0 to + 107 and the 
remaining from 0 to + 57. 
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8. Dii^ital output 

12 sets (16 bits each) of digital output terminals. 

Two sets are of electronic contact operate (EGO) type 
and the rest are of voltage sense type, 

9- Process interrupt switches : 5 sets of 16 bits each of 
contact sense type. 

10. Interval Timers ; Three interval timers, A,B and C. 

Timer A and B are wired for .25 ns and 1 ns time bases 
respectively. Timer C acts as a z-eal time clock, 

A-2 AQ 20 Specifications fl] 

This is an ICIL make analog computer, 

1. Total 20 operational amplifiers (+ lOV) with five inputs 
(3 of unity gain and two with a gain of 10). All of them 
can be used as summers /inverters while 10 of them can be 
used as integrators also. 

2« 4 Quarter square multipliers (+ 10 V) 

3. One universal diode function genei’ator (+ 107). 

4. 30 ten turn hand - set helical potentiometers. 
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C this SU8R0U~ir-.b PETS Tf-E ANALOG CCPPUTE 

C IL RESET YCDE 

data JDATA(l),„':ATA(2),vtATA(3>/6,127,2/ 
call CD ( 11 00 1 , ..iJ at a ( 1 ) , V C a T a ( 3 ) ) 

RETURN 

ENC 

a ST ORE KcSeT 

ar.CENC 


/ / J C E 
/ / DUE 

flEELE’'E SbJ 

// FOR SST 

subroutine SST UTST ) 

CCKHON BALU,KF.KS 

C JTST WILL BE 1 IF CONFUTE TIME IS OVER 

^tst=balu 

RETURN 

Ei\r: 

SST 


«STORb 

«CCE!\i; 


LU > 



// vj"? 

/ / D I ? 

«IEl.E-E ■ M = 

// FCE ^JrE 

Si.:" RCU" T.E 
CCF’-D?'.' 

c Tf-rg PQir 'r, ■■■ E E:‘=;tFEE AFThF 

C «• r E b L I! ■ S 3 E A : - 1 ^ i-' E N C C i“ F L T F 
KF = kP4.i 

IF(KP-KS) 10, ^'^>12 
■C CALL TILS" 

PEC S„i R Sx 

1? BALL=1 

RET URL 
E?.. D 

-STOP" ; .vj- 

«ccE‘ i: 


// JCF 
// DiJp 

«CElFTp : ; 

// FCK T'F'EF 

SlRPOU'!‘!LE TTN.^r 


C THIS IS I'D Pui 

THE analog computer ? ■. CCMFI.TE 

C * ® M C n E FOR NCR: 

:HAN SECONDS. 

external time 


CALL T!MER(TI“= 

. k f 25 L C C ) 

RETURN 


END 


«STCF?t ' 

! ME F 

fiCCENC 



EVERY ,;-T SECC 

T f y ^ ^ Q £ C 


n jc& 

// DUP 
// FCR IDLE 

subroutine IDLEUIf^) 

C 6* THIS ROUTINE IS iNTENEEC TC PROVIDE aNY DELAY REQUIRED AT ANY 
C DELAYS AREREQUIKEC, FCR EX AnFlE> AF TER A hCCE CHANGE CQHNAND TO 

C COMPUTE SINCE NC CQNFLTATICNS SHOUlD BE CGNE UNTIL ANALOG HASC 

C ®« MCDF AMD IT TAKES AbCLT iCCMSfcC FCR MECH CONTACTS TO CLOSE 
C THE LOOP BELOW IS A £.5MS LOOP 

N=TIM/,0085 
CC1!=1,M 

G2=1.0/l,®2./2./l./3./2. 

1 continue 

RETURN 

END 

^DELETE IDLE 

«STCRE lULb 

«CCENC 



' H i S i S I ■ ~ Y 't; ^ ” C F ^ C G P A X 


h R R i G S > 3 - L G L 

'Em ;- r, 

P P V 3 p 

( Y P ~ k’ R ■ P P » < E V ri _ /t K L } 

R E A ^ PI j K 1 » G 3 '- > K 1 » K F » M 1 » M 2 

C I '•■'E MS I CM X ( • Y1 ( c.^ ) 3 X2 ( E5 ) 3 Y3 3 5 

Ynl=10 2. 5 

P I = E . 1 4 1 (5 

•VSj = PI / 4 . ■”• ( ( 0 . :js2 A-i-;2 ) «S2 . «113 . + ( C 
VS ? = P I /4 . -s- { ' i; , i. 26 3-y « 2 ) « IS , « 2i 5 . + ( 

1 + C G .0355**2 ) •"' 2 6 . * k 2 a . 'll 
K 1 ~ 5 6 , 5 
K2--^'5GG'D . 

HKCS = 626 . 7 
HCS=247 . 64 
-=n . 00235 
P 1 = 0 . 8 

C= 0 . 0891 *SQR ■ 1 U . 3C 4>; ) 

CF5 = 600n./(l.J7.i34*-sC.6*2C . ) 

P S T = 6 6 ■< . 5 9 5 *Ji 
RCST1=76 .5111 
TST1=36Q. 

TH3=1.Q73033 
TR1=1.3P4795 
TH2 = 1 . 76 
1P2*1.3795 
HS‘^? = 820. 
hCST 2 = 40 .2-576 
TST2=54n . 

KF=tOR. 5/110030 . 

EAV=0. 0392 
VEL=3750.»0,C2UU4/(9.*FI«CAV**2»224.? 

TCL=200./VEL 

CS1=0.164 

CS2=0.164 

M1=F 1*0 . 0523*4 . 0. an8«32 .•all3 . -sSOFO . 

1 + F 1*0 . 04*2 . -so . '-'6 45 ^112 .*22 6 . *8020 . 

M2 = P 1*0 .04*2 . «U . 04^*C .1*3C .*440 . «802Q PI :.C . 0392*7 . -"-O . 009*26 
1*2 24. *80 20, 

ChS=:18660.0/(ll4.2**C.6*2C. ) 

G SHI = 60 00 , 

QSM2 = 18660 . 

RCS1=77.4327 
HS1=663 . 786 
hS2=8l9. 112 
RCS2=4a.l474 
WC1=107.334 
TSH1=380.099 
TSK2=568.09 
A1=0 .3041 
A=A1 





j. , ! 

Q2 =6 55 .756 
G. 3 = -2 6. 147 4 
G4=F03 12 

G5=il3.334 
Gfi=376.799 
■37=5 58 .09 
H1=0. 6 
r2 = & ,8 
= n .4 
h- ^ = 2 , 0 
H5 = l .0 
H 6 = fj , 4 
H7 = :' .0 

X(3 )=X(1>-G1 
X (? )=X( 2) -G2 
X(:5)=X<3 )-Q3 
X(4)=X(4)~G4 
X(5)=XC5)-G5 
X(6) =X(6 )-G6 
X ( > ) = X ( 7 ) -G 7 
in KBG=KBG+1 

XCl)=X(l)+Gi 
X(2)=X(2)+G2 
X(3)=X(3)+G3 
X( 4) =X (4 )t-G4 
X(5)=X(5)+G5 
X(6) =X(6 )+G6 
X(?)=X(7)+G7 

Xl( I )=(HR1-X(5)+WDS1)/VS1 
PS1=(K1«X(?)-K2)«X(1) 

TSl=THlft(X(2)-HbTl)+TPl«(X(l)-RCSTl ) + TSTt 
QS1=CPS*X(5 (X (6 )-TS1) 

F11 = (QS1 + HR1»HWQS-X (5 )«X (2 ) + WDSl*HES) /( X ( 1 )«VS1 ) 
X2( I) = ( F11+ J^^Pbl«6. «X 1{ I ) / ( X( 1) «X (1 ) 5 ) / (1 .C-K 1-^ 
PS2=(K1«X(4)-K2}«XC3) 



1 = C * ( p s 2 i? /. •' ' •• ) -«■ . : 

>' 3 «' : ) = ^ X ; 5 ) " ;:s 2 ) / s 2 

’S2--'H2';f-(X {/- )- -bl> (X (2 )-P'.S''? 5 -^5^2 

■p S 2 = C H S * A 1 ■«' « [•. 1 '"■ ( X ( y ) - ' S 2 ) 

2 2 3 = ( Q S < + y ' 5 ) : 2 ) 2 « !- 2 S - A 1 •”• X ( 4 ; ) / , / ( 3 } «- v 

X ^ ( : ) = ( ' 3 3 + >; «■ P c > 6 . ■■■ X 3 ( I ) / ( X < 3 '/ * ( 2 ; " } 

X 5 r : ) =i/ P/ - Cl « r r i ? - F 3 2 ) - X ( 5 ) / T CL 
X 6 f I ) = i Q3 1 -• C S : ^ / C •- 5 1 « M ) 

X7 ( : ) = f CSX2-CS2 ; / ( Cb2«^'2 ) 

X (2 ^ = -y 2{ r> /H 2 
x(l)=-X3{ I )/H5 
X ( 4 ~ X /, ( T) / ■•-! 4 

X(5)=-X5f I 
tit') =-X6( ; )/. 6 
X ( 7 ) =-X7( ■ ) /H / 

3 C 1 nil t Al = 2 » / 
in:' X ( K /■ L ) = X ( K A L : / 4 . -i 

V.F ' =0 

i: c ‘^■Q K = 1 , 6 

K R = K + 1 
V ALIE=K (KR) 

CALL k R A C ( V A L. t f ^ p ! ) 

4C .,-FT=JP'^ + l 

IF(KL) 21,21,2b 
21 CALL CO^^PinCL.! 

KL = KL + " 

CALL IDLE(3.G) 

K N = 4 1 1 2 

25 CALI. RDAKX.KL ) 

X(1)=X(1)«H1 

X f2 ) = X< ?)■»'-’? 

X(3)=X(3)*H3 
X{4)=X(4)»H4 
X(5) = X(5)<»H5 
X ( 6 ) =X ( 6 )*H6 
X(/)=X(7)«H7 
XI { i ) =X (1 ) 

X2 ( 1 )=X(2) 

X3( I )=X(3) 

X4( I )5X(4) 

XSC I )=X(5) 

X£'< I ) =X (6 ) 

X7 ( I ) = X{ 7) 

IFCKBG-KBG/lD^y f ) 12,13,12 
13 1=1+1 

12 CALL SST(JTST) 

IF(I-IOO) 26, 2b, 28 

26 IF(JTST) 10,10,28 

2g CALL RESET 

CC 60 KS=1/ I 
XI (KS) =X1( KS )+Gl 
X2(KS)=X2£KS)+G2 
X3(KS)=X3(KS) + G.< 

X4 fKS) =X4£ KS ) + G^ 

X5(KS ) = X5(KS)+Gtj 


m 



m \ji m C/1 c/i 


) = A t i 3 ) 1’ t 
} = X 7 ( . 


c. , ' 0.: ) X - f KC ) , ( KC ) » X3 ( KC ; . x"-; Kc ; 
"C^'‘'AT(4(2/, 

:CM IfAJE 
: A M V ! A n 


G ■! , X 7 ( s<r ( 


!r( -’fER hYSR 


sCCE 

sSFCCr RAGS HYBR 
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